NACA TN 3567 


govt. doc. 


y?..H 2i/r: 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3567 business and 

technical dept- 


r 


STUDY OF SCREECHING COMBUSTION IN A 6 -INCH 
SIMULATED AFTERBURNER 

By Perry L. Blackshear, Warren D, Rayle, and Leonard K. Tower 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

October 1955 


1 


3070 


NATIONAL ADVISORY COMMITTE FOR AERONAUTICS 


TECHNICAL NOTE 3567 

STUDY OF SCREECHING COMBUSTION IN A 6-INCH SIMULATED AFTERBURNER 
By Perry L. Blackshear, Warren D. Rayle, and Leonard K. Tower 


SUMMARY 

As part of a general research program on screeching combustion at 
the NACA Lewis laboratory, an investigation was conducted to develop 
screech instrumentation and to study the mechanism of screech in a 6- 
inch-diameter simulated afterburner. 

Probe microphones were developed that can be utilized to measure the 
frequency, relative amplitude, and relative phasing of the pressure os- 
cillations at various positions within a screeching combustor. In cali- 
brating these microphones to determine absolute values of pressure am- 
plitude, a new theory was proposed to account for the nonlinear attenua- 
tion of high-amplitude sound in tubes. 

The acoustic oscillations accompanying screech in the 6-lnch after- 
burner consisted of the first transverse (sloshing) mode in the hot gases 
downstream of the flameholder. 


INTRODUCTION 

The trend toward higher flight speeds for military aircraft creates 
the need for jet engines of Increased thrust. The trend in afterburner 
and ram-jet combustor development has therefore been toward combustors 
that operate at higher pressures, higher combustion temperatures, and 
higher velocities. With the development of these high-output combustors, 
the phenomenon known as screech has been widely encountered throughout 
the aircraft industry. Screech derives its name from the high-pitched 
audible sound accompanying the phenomenon. Other manifestations of 
screech are high-frequency pressure oscillation in the combustor and an 
increased rate of heat transfer, which have resulted in rapid deteriora- 
tion or failure of combustor shell, flameholder, and other combustor 
parts. ^ In most Instances, therefore, screech is undesirable and is to 
be avoided. It is essential that an understanding of the screech phe- 
nomenon be obtained and that a method of controlling or eliminating 
screech he devised. 
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Screech has received some detailed study in rocket (refs. 1 and 2), 
simulated ram- jet (ref. 3), and afterburning (ref. 4) engines. In ref- 
erences 1, 2, and 4, screech was identified as a transverse mode of os- 
cillation hy con 5 )aring the experimentally observed frequency with the 
frequencies calculated from the various acoustic modes possible in an 
idealized chamber resembling the combustors employed. In reference Z, 
screech was Identified as a transverse mode by schlieren pictures of the 
motion of the flame near the flameholder and by inspection of the fre- 
quency obtained. A more detailed study of the screech mechanism is 
desirable, since (l) a wide assortment of idealized modes can give evi- 
dence similar to that reported in references 1 to 4, and (2) conditions 
within the combustors suggest a much more complicated picture of the 
oscillation than the idealized acoustic treatment admits (ref. 5) . 

The Investigation reported herein, which is part of a general pro- 
gram to investigate screech at the NACA Lewis laboratory, hah. a two-fold 
objective: (l) to develop suitable instrumentation for studying screech' 

in an operating combustor, and (2) to make the measurements of screech 
necessary to assess the screech mechanism. The investigation was con- 
ducted with a 6 -inch-diameter direct -connect combustor. The combustor 
configuration and operating conditions simulated a turbojet afterburner. 
The instrumentation, which was developed and evaluated for studying 
screech. Included a probe microphone capable of measuring the pressure 
fluctuations at any point within the combustor. This permitted a detailed 
study of the spatial distribution of amplitudes, frequencies, and rela- 
tive phases. Ionization gaps were utilized to indicate the position of 
the flame front in the vicinity of the flameholder. In order to calibrate 
the probe microphone, an experimental and theoretical investigation was 
made of the attenuation of large-amplitude sound in tubes. 

From surveys of the screech-amplitude distribution, frequency, and 
relative phasing, the oscillatory motion of the gases in the combustor 
was deduced by comparing the experimental results with phemomena pre- 
dicted from theoretical consideration of various modes of acoustic os- 
cillation. A limited study was made of the effect on screech amplitude 
ana frequency of such system variables as combustor pressure, air-flow 
rate, fuel-air ratio, and flameholder type and position. 

The work reported herein was conducted during 1952-53, and has been 
reported with additional information in an NACA classified document. 


APPARATUS ANT) PROCEDURE 
Installation 

A dlsigram of the combustor Installation is shown in figure 1. The 
inlet and outlet ducts were connected to the laboratory air supply and 
exhaust systems. Air -flow rates and combustor pressures were regulated 
by a remote-controlled valve upstream of the combustor and by a variable- 
area combustor exhaust nozzle. A tubular production-model turbojet 
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combustor was installed in the duct upstream of the research comhustor 
to prov^ide conditions in the inlet stream simulating conditions at the 
inlet of a turbojet af terbtirner . 


Research Combustor 


The research combustor (fig. 2) had an inside diameter of 6 inches 
and a length (measured from the position of the flameholaer to the posi- 
tion of the exhaust nozzle) that could be varied from 34 to 45^ inches 
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by •brav'ersing tbe riameholder. The riameiiolder was -tx-aversed axially by 
means of a vire cable and a winch. The combustor having the variable - 
area exhaust nozzle and the water-jacketed combustor shell shown in fig- 
ure 2 was used for only a part of the investigation reported herein; 
during the early phases of this investigation, the combustor had a fixed- 
area exhaust nozzle and an air-cooled outer shell. Fuel, JP-4 or gaso- 
llnG^ was supplied to the combustor through an air— atomlz ing spray bar 
(fig. 3) located 100 Inches upstream of the exhaust nozzle at the posi- 
tion indicated in figure 1. Four flameholders were investigated: two 

single diametric V-type flameholders 1^ ana 2j inches wide and two 60°- 
single-cone flameholders with area blockages of 30 and 46 percent. 


General Instrumentation 

Air flow was metered by a concentric-hole, sharp-edged orifice in- 
stalled according to A.S.M.E. specifications. Fuel flow was metered by 
calibrated rotameter . Combustor— inlet pressure was measured by four 
wall static taps as indicated in figure 1. Combustor-inlet ten^jerature 
was measured by chromel-alumel thermocouples connected to a self- 
balancing recording potentiometer. A thrust target was en^Jloyed to 
measure the momentum of the exhaust. This was later replaced by a water 
quench spray, followed by a thermocouple array. The rise in enthalpy 
across the system then indicated the thermal efficiency. 


Special Instrumentation 

Microphone probe development . - An instrument was required that 
could be used to measure the frequency and amplitude of the pressure 
oscillations acconq)anying screech and to determine the pressure distribu- 
tion of any standing wave system in the combustor. To determine the 
pressure distribution, two basic techniques were apparent: (l) move the 

sensing element to various positions in the combustor, or (2) utilize a 
movable probe to pick up the pressure signal at various locations and 
transmit 'this signal to the sensing element located outside the combustor. 
There exist water-cooled pressure pickups that can withstand limited ex- 
posure to flame when mounted on the conkiustor wall. The life of these 
pickups is in constant jeopardy with designs einployed so far, and to 

adapt these pickups to immersion in the hot gas stream would hardly seem 
feasible. 
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Therefore^ probe microphones for measuring the freq.uency and rela- 
tive amplitude of screech in an operating combustor were designed^ fab- 
ricated^ and subjected to preliminary evaluation by measuring screech in 
the 6-inch-dlameter simulated afterburner. The instrimients showing the 
greatest promise were then selected for more detailed calibration and 
additional studies of screech. The necessary calibration consisted of a 
measurement of the attenuation of pressure waves between the open end of 
the probe and the microphone. 

Description of instruments . - The probe microphone is shown in fig- 
ure 4. A water-cooled probe that could be traversed diametrically through 
the combustor was connected to a 50-foot coil of l/4-inch copper tubing 
that behaved acoustically as an infinite tube in that it produced neg- 
ligible reflection of pressure pulses entering the tube. The sensing 
element^ a microphone^ was connected to this "infinite” tube by a short 
l/32-inch-diameter lead at a point near the combustor. To prevent com- 
bustion gases containing soot and moisture from entering the probe during 
normal operation^ a small flow of clean^ dry air was maintained through 
the "infinite” tube into the combustor. The air purge was not adequate 
to keep the tube completely clean^ since rapid changes in combustor 
operating conditions caused pressure surges that sometimes carried com- 
bustion gases into the tube. To avoid changes in the tube attenuation 
constants resulting from soot deposited from these combustion gases ^ it 
was deemed necessary to use a tube of reasonably large diameter ( 3/16- 
in. I.D.). 

Special probes were fabricated to check the reliability of the phase 
measurements obtained with the standard probe microphone. These "phase" 
probes were intended primarily to produce reliable phase measurements; 
the reliability of their amplitude readings was assumed to be small. A 
diagram of such a probe is sho^wn in figure 5. The distance from the end 
of the probe to the microphone was minimized to avoid variations in in- 
ternal phase lag. The bore of the tube was kept small to enable viscous 
damping to minimize the effect of reflected waves. An ample supply of 
cooling water was passed through the jacket of the probes^ keeping the 
gas within the small tube at a temperature above the dewpoint of the 
combustion products . 

Readings obtained with the phase probes confirmed the validity of 
the phase measurements made with the standard probe microphone. The 
phase measurements are therefore presented under RESULTS AND DISCUSSION 
vmthout qualification as to the instrument employed. 

Calibration . - In order to calibrate the special screech instruments 
for frequency and amplitude response^ the apparatus shown in figure 6 was 
utilized. The apparatus for generating low-amplitude pressure waves 
(fig. 6(a)) consisted of a resonance chamber driven by a piston speaker. 
The apparatus for high- amplitude pressure waves (fig. 6(b)) consisted of 
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a similar resonance chamber driven by a high-pressure air Jet that was 
admitted intermittently by a rotating "chopper" disk. The readings of 
the research Instruments were compared with the readings of a reference 
microphone that had previously been calibrated with a known steady- state 
variation in pressure. From these data the attenuation coefficient for 
the probe microphone was computed. 

A theoretical analysis was made of the damping of high -amplitude 
pressure waves^ and the resvilting theoretical equations were utilized to 
correlate the experimental data on attenuation of high- amplitude waves. 
Theoretical considerations were also utilized to correct for the change 
in amplitude accompanying the transmission of a pressure wave across a 
temperatiire interface such as that existing between the hot gases in the 
combustor and the cold purge air flowing through the probe microphone. 


Attenuation of Low-Amplitude Sound in Tubes 

The attenuation and distortion of periodic pressure waves in tubes 
must be considered in designing and calibrating the probe microphone. 

The apparatus shown in figure 6(a) was utilized to measure attenuation 
constants for low-amplitude pressure waves and for various probe lengths. 
The so-called Kirchoff attenuation constant for l/4-inch tubes was 

computed from the relation 


^2 

X2 - xp ( 1 ) 

where the symbols ^ve the mearri.ngs listed in appendix A. The experi- 
mental values obtained from equation (l) were approximately 15 percent 

greater than the theoretical Kirchoff value; the data are presented in 
figiu'e 7. 


Attenuation of High-Amplitude Sound in Tubes 

Screech oscillations having a peak-to-peak pressure amplitude as 
great as 20 Inches of mercury have been encountered. For such pressure 
oscillations the magnitude of the attenuation coefficient is reported to 
be dependent on the amplitude (refs. 6 to 8) . Since these references do 
not agree on the shape of the attenuation-against-amplitude curve, the 
data contained in the references were not considered adequate for com- 
puting the damping of high-frequency, high-amplitude waves in the l/4-lnch 
^ experimental and theoretical investigation of the damping of 
high-frequency, high-amplitude pressure waves in small tubes was there- 
fore undertaken as part of this Investigation. 
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Theoretical . - A derivation of equations describing various factors 
affecting the behavior of high-amplitude sound in tubes is given in 
appendix B. The factors can be briefly summarized as follo\/s: 

The sound pressirre in the entrance of the probe is related to that 
of the soiirce by the equation 


Pr - Pe = H Pe 

It is well-known (ref. 9) that high- amplitude sound waves propagat- 
ing in free air tend to develop a saw-tooth wave form; that is , the \mves 
generate higher harmonics. The amplitude below which the n^^ harmonic 
of a wave of frequency f will not grow is 


~ 0.29'^ can"^(n-l) 
Pmax “ f(n+l) 


(3) 


where 

0.2921 depends on y .= 1.4 

a damping constant for fundamental component 

m exponent of f in a (a = a’f^) 

In free air m = 2 and a is rather small for the frequencies of 
primary interest. For viscous damping in tuhes^ m is l/2; and^ if a 
is defined as Vf. equation (3) becomes 


^max 


0.292 ca’ 


In (n-1) 

'Vf tth 


(4) 


For the case where this value is exceeded for the first several harmonics, 
the length of propagation for a sound wave to become saw-toothed is 


X 


~ 


= In 



e .max 




max 


- 0 . 


292Xo^ J 


(5) 


If the amplitude-dependent attenuation coefficient is defined as 


dp 

max 
dx _ 

Pmax 


( 6 ) 
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this constant can be expressed as 


(again for y = 1.4). 

A plot of equation (7) is presented in figure 8. Values of the 
amplitude-dependent damping coefficient Op were combined with the 
experimental values of the Kirchoff damping coefficient and the 

resulting expression for the combined damping effects was integrated to 
obtain the curves shown in figure 9, where values of p^^^ axe plotted 

as a function of tube length x for three frequencies. Figure 9 indi- 
cates that even for high amplitudes the Kirchoff damping has a rather 
large effect in determining the over-all attenuation. 

Therefore, the factors affecting a high-amplitude wave as it enters 
and propagates in the tube are as follows: First the amplitude is re- 
duced on entering the tube according to eqiiation (2). If the wave form 
is sinusoidal, it will become saw-toothed after propagating a distance 
given by equation (5), dioring which time it is attenuated by Kirchoff 
damping (eq. (l)). Once the wave becomes saw-toothed, it is attenmted 
according to the integrated eqioation shown in figure 9. Consideration 
of these factors yields the theoretical curves shown in figure 10. 

Experimental . - Figure 11 shows the wave form recorded at 1080 
cycles per second for various tube lengths in the probe microphone using 
the apparatus of figure 6(b) . From these wave forms and from simulta- 
neous voltmeter readings, the mean-to-peak amplitudes were deteimined. 

Figure 12 presents a comparison of the pressure amplitude indicated 
by the probe microphone with that indicated by the calibrated reference 
microphone for frequencies of 1080, 2400, and 3000 cycles per second. 

Also presented in figure 12 are the theoretical values of figure 10. 
Values of Op obtained from the curves of figure 12 are compared in 
figxire 13 with the theoretical values from figure 8. Also shown in fig- 
ure 13 are values of Op obtained from references 6 and 8. 

Discussion , j- From the foregoing results it is seen that probe 
microphones used in measuring high-amplitude, high-frequency sound become 
less sensitive as the probe length is increased and as the wavelength 
decreases. Fortunately, observed screech frequencies tend to decrease 
as combustor diameter is increased (e.g., ref. 4). The resulting design 
principle is fairly obvious; for maximum sensitivity of the probe micro- 
phone, the distance from the end of the probe to the microphone should 
be kept small. It is possible to obtain a reasonably accurate picture 
of the wave form by the use of a probe microphone only if the wave 
travels less than a distance x* as defined in equation (5). 


+ Pmax)^ - 1 


(7) 
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The agreement between the observed and computed nonlinear attenua- 
tion coefficient is surprisingly good. If calibration facilities are 
not available, these values may be used for computing the amplitude- 
dependent attenuation until some more refined method becomes available. 

References 6 to 8 postulate turbulence as the cause of the nonlinear 
attenuation of high -amplitude sound in tubes. Figure 13 shows that the 

theoretical nonlinear damping coefficients derived in appendix B are in 

reasonable agreement with both the experimental values of references 6 
to 8, which were obtained in large tubes, and the experimental values 
obtained in this investigation with l/4-inch tubes. Consequently, it 
appears that Op is independent of tube size and hence is not a direct 

function of Reynolds number. The theory presented in appendix B is 

therefore believed to afford a better explanation of the observed non- 
linear attenuation. 


Effect of Temperature Discontinuity on Probe-Microphone Reading 

All the foregoing effects have been for an assumed constant mean 
temperature. When the gas temperature within the probe differs from the 
temperature in the combustor, a reflection will take place where the 
sound wave enters the tube. Assimiing that particle displacement and 
sound pressure are eqiial on either side of this interface, the ratio of 
the transmitted wave to the incident wave amplitudes can be computed as 
a function of the temperature ratio of the gas on either ^icie of this 
interface. This ratio is given in the following equation: 

( 8 ) 

Pine 2 . + 

1 ^inc 


Flame Detection 

Ionization gaps were used to detect the local motions of the flame 
during screech and to give an indication of particle motion in a com- 
bustor. The ionization gaps used are shown in figure 14. The use of 
these gaps coupled with a surface -mounted pressure pickup provided an 
instantaneous picture of the relation between the local pressure and the 
flame-front position near a flameholder. 

Several variations of photocell probes were made and tried. The 
results were not satisfactory. A periodic variation in the limiinosity 
of the flame zone downstream of the flameholder was obtained with the 
photocell probe sketched in figure 15. This probe consisted of a small 
periscope purged with clean, dry air flowing through the periscope into 
the combustor. The photocell recorded an integrated value of the light 
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emitted within a narrow cone between the periscope and the wall. The 
photocell probe produced a periodic signal during screech and therefore 
shows promise; it has not been used sufficiently^ however^ to warrant 
its recommendation. The ionization gaps^ on the other hand^ which were 
simpler to construct and to operate^ have given very interesting data 
and are recommended for further screech investigations. 


Procedure 

The relative screech amplitude at various positions across the 
diameter of the 6 -inch afterburner was measured at different distances 
downstream of the flameholder. The distance between the flameholder and 
the probe microphone was varied by changing the axial position of the 
flameholder. The pressure amplitude at the wall of the combustor was 
measured at 18 points along the length of the combustor by leaving the 
flameholder fixed and moving the microphone to the various positions of 
measurement. During the process of making these sui^eys of relative 
screech amplitude^ the frequency spectrum of the sound was also recorded 
at various positions in the combustor. 

The relative phasing of the pressure signal at different positions 
in the combustor was determined with two microphone pickups. Phase 
measurements were made between a reference pressure signal and the pres- 
sure signal obtained in vertical diametric surveys at three distances 
downstream of the flameholder 4^ and 7-i in.). Phase relations were 

also measured at 18 points along the length of the combustor. With two 
ionization gaps located 180^ apant^ 3/4 inch downstream of the flame- 
holder^ and extending 3/4 inch from the burner walls ^ oscillations of the 
flame were detected. Simultaneous measurements with the ionization gaps 
and a microphone pickup were used to determine the relative phasing be- 
tween the flame position and the pressure in the combustor. 

The peak-to-peak pressure amplitude was measured at various pressures 
and air-flow rates in the 6-inch af terbirrner . At each of several air-flow 
rates the combustor pressure was progressively lowered to determine the 
limiting pressure below which screech was not encountered. Measurements 
of screech amplitude and frequency were made at various fuel-air ratios 
in the combustor; these data were obtained with the fixed-area exhaust 
nozzle^ however^ and changes in fuel-air ratio were accompanied by var- 
iations in combustor- inlet velocity and canbustor pressure. 
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RESULTS AND DISCUSSION 
Screech Mechanism 

Amplitude distribution . - The relative screech amplitude at various 
positions across the combustor diameter is shown in figure 16 for sur- 
veys made at different distances downstream of the flameholder. The 
data of figure 16 were obtained by traversing the flameholder axially 
along the length of the combustor. At certain locations of the flame- 
holder in the combustor, the longitudinal mode of oscillation would 
sometimes appear, and the indicated unfiltered amplitude of screech 
would rise. Figure 17 shows am oscilloscope trace of a combined Intense 
low-frequency note with screech and a phase relation of the two frequen- 
cies obtained with appropriate band- pass filters. There Is an apparent 
modulation of the screech amplitude by the low-frequency note. For most 
positions of the flameholder, the 3300-cycle screech constituted the 
dominant oscillation. For other positions, however, the low-frequency 
490-cycle oscillation was dominant, as indicated in figure 16. Figure 
18 shows typical spectra of sound recorded at the combustor wall for both 
the low- and high-frequency oscillations. 

In all surveys shown in figure 16 where the 3300-cycle oscillation 
was dominant, the amplitude exhibits a pronounced dip near the center of 
the duct. This is the amplitude pattern to be expected with a transverse 
mode of oscillation, as shown in figure 19. Such a pattern could not 
result from radial or longitudinal oscillations. 

Figure 20 shows a comparison of a sound spectimm at the combustor 
wall with a spectrum at the center of the duct for conditions producing 
the amplitude dip at the center. Only the high-frequency part of the 
spectrum (3300 cps) is attenuated at the center. It is concluded, there- 
fore, that the frequency associated with the transverse mode is about 
3300 cycles per second. The data shown in figures 16, 17, 18, and 20 
are raw data taken with 3 feet of transmission line piping the signal to 
the microphone and subject to the corrections previously discussed. 

For two positions of the flameholder, the amplitude at the wall was 
measured at 18 points along the length of the combustor with a single 
microphone moved from point to point to obtain the data shown in figure 
21. These data show a peak in the screech amplitude about 5 inches 
downstream of the flameholder. This peak lies in the region where the 
TnflxiTmiTn heat-release rate might be expected. 

The circumferential distribution of sound at a plane 1 ^ inches 
downstream of a V-gutter flameholder is shown in figure 22. These data 
were gathered in three successive tests, the flameholder being rotated 
45° between tests. The filtered screech signal is plotted against posi- 
tion relative to the gutter. A strong maximum is apparent at locations 
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perpendicular to the gutter^ thus indicating that the gutter orientation 
controls the direction of the transverse oscillation. Also plotted is 
the higher-frequency content of the signal; this, too, is minimum behind 
the gutter but does not demonstrate so marked a peak as the screech 
f vindamental . This may be interpreted as suggesting a nonsinusoidal driv- 
ing force which acts over a relatively large portion of the cross section. 

The data presented in figures 16, 21, and 22 are considered con- 
clusive evidence that the screech oscillation studied comprised a stand- 
ing transverse wave. The circumferential and radial distribution follows 
the general shape demanded by such a wave; the longitudinal survey may 
also be so interpreted. The amplitude peak shown on figure 21 does not 
indicate the presence of a standing longitudinal wave; rather it signi- 
fies the location from which the transverse oscillation is driven. If 
a combined longitudinal and transverse mode existed, the frequency would 
be 



( 9 ) 


where njl is the reciprocal of the half wavelength of the longitudinal 
oscillations, and d is the burner diameter (ref. lO) . If this longi- 
tudinal component is to have an appreciable effect on the indicated fre- 
quencies, the half wavelength would have to be of the order of a diam- 
eter. With a half wavelength of this magnitude, multiple pressure nodes 
would occur in the 38-inch length surveyed. Since multiple nodes are 
not evidenced in figure 21, there exist no longitudinal components which 
can appreciably affect the frequency. The screech frequency can there- 
fore be expressed by the relation 


f 


c3 

d 


( 10 ) 


The values of p (ref. 10) for the first two pure transverse modes of 
oscillation are 0.568 and 0.972. 

Frequency . - The observed screech frequency of approximately 3300 
cycles per second can be produced by the first pure transverse mode of 
oscillation (P = 0.568) if the speed of sound is 2920 feet per second, 
or by the second transverse mode (p = 0.972) if the speed of sound is 
1710 feet per second. These values for the speed of sound correspond to 
average temperatirres of 3850° and 1250° K, respectively. The higher of 
these temperatures approaches the theoretical flame temperature; the 
lower approaches the temperature of the inlet air. Consideration of 
screech frequency therefore indicates that the oscillation could consist 
of the first transverse mode in the hot gas or the second transverse 
mode in the inlet stream. For transverse modes of oscillation, the 
pressure will be out of phase on the opposite sides of the combustor for 
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odd modes and in phase on opposite sides of the combustor for even modes 
of oscillation^ along any diameter except an antinodal diameter. Phase 
measurements should therefore make it possible to determine whether the 
first or the second transverse mode of oscillation occurs. 

Phase measurements . - By means of the special phase probes^ vertical 
diametric surveys of pressure amplitude were made at three distances 

•z n 

downstream of the flameholder {-g, and 7^ in.) to obtain the data 

shown in figures 23 and 24. ’A typical oscilloscope trace is shown in 
figure 23. The phase of the pressure signal relative to the signal ob- 
tained 3/4 inch downstream of the flameholder at the top wall of the 
combustor is indicated in figure 24. The 180° phase shift across a 
diameter indicates that the first transverse mode rather than the second 
mode of oscillation persisted in the combustor. The data of figirres 23 
and 24 Indicate that there is also a phase shift along the length of the 
combustor. The relative phase of the pressirre signal along the top of 
the combustor is shown in figure 21. These data also indicate the ex- 
istence of a pronounced phase shift in the vicinity of the flameholder^ 
signifying that in this region the waves have longitudinal components. 

The absence of a phase shift in the hot gases downstream of the flame- 
holder indicates the existence of pure transverse oscillation in this 
region. 

Circumferential phase relations were sought. The results were not 
conclusive and could be interpreted as evidence either of standing or 
traveling transverse oscillation according to the desire of the inter- 
preter. This may be attributed to the large separation (90°) of the 
instrumentation positions, to the nonsinusoidal nature of the wave, or 
to a phase lag between the driving force and the oscillation. 

With two ionization gaps located 180° apart, 3/4 inch downstream of 
the flameholder and extending 3/4 inch from the burner walls, the signals 
shown in figure 25(a) were obtained. The flame striking the ionization 
gaps was 180° out of phase on opposite sides of the combustor. At posi- 
tions farther downstream a random component was superimposed on the 
periodic component, making the signal difficult to interpret. 

Simultaneous measurements of pressure and flame -produced ionization 
at a station 3/4 inch downstream of the flameholder (fig. 25(b)) showed 
that the outward flame-front motion is in phase with the pressure at this 
station. On the basis of these phase measurements, a time sequence of 
events in the screech cycle can be constructed as shown in figure 26. 

The relation between the pressure and the flame-front displacement is 
precisely that which wo\ald accompany a transverse oscillation. 
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Effect of System Variatles on Screech 

Pressure . - The effect of pressure on screech amplitude is shown in 
figure 27. The peak-to-peak pressure amplitude showed a general Increase 
with increase in combustor pressure. These measurements were made with 
a flush-mounted microphone located upstream of the flameholder and there- 
fore do not represent the maximum values of amplitude (see fig. 2l) . The 
air-flow rate is noted for each data point in figure 27. At an air-flow 
rate of 1.5 pounds per second, the pressure limit for screech was 20 
inches of mercury. Detailed measurements were not made of variations in 
frequency with pressure; sufficient observation was made, however, to 
ascertain that no marked change in frequency occiirred as pressirre was 
varied. 

Through-flow . - The combustor screeched throughout most of its 
operating range, as shown in figure 28. Screech was encountered at pro- 
gressively lower pressures as through-flow was decreased. During screech 
the combustion efficiency was 90 to 100 percent; the efficiency dropped 
abruptly by approximately 35 percent when screech ceased. 

The limited data of figure 27 show no consistent variation of 
screech amplitude with variation in flow rate. No effect of flow rate 
on screech frequency was noted for the limited range investigated; these 
data are presented in table I. These tests were conducted with a fixed- 
area exhaust nozzle, however, and the maximum variation in ccmbustor- 
Inlet velocity was about 20 percent. 

Fuel-air ratio . - The effect on screech of variations in fuel-air 
ratio is shown by the data of table I. The fixed-area nozzle was used 
to obtain these data; this did not permit changes in fuel-air ratio with- 
out simultaneous variations in combustor-inlet velocity and combustor 
pressiire. Only qualitative trends can therefore be noted. Fuel-air 
ratio exhibited no pronoimced effect on screech amplitude. The effect of 
fuel-air ratio on screech frequency was as follows. Maximum frequency of 
3500 cycles per second was obtained at an equivalence ratio of 1.146. 
Screech frequency decreased as fuel-air ratio was decreased below this 
value. The frequency was 3350 cycles per second at an equivalence ratio 
of 0.934 and 2800 cycles per second near the lean blow-out limit. When 
the fuel throttle valve was abruptly closed, combustion persisted for one 
to two seconds . During this time the screech frequency continuously 
declined to a value much lower than 2800 cycles per second. 

Flameholder type . - Two single diametric V-type flameholders and two 
60° single-cone flameholders were investigated. With the V-type flame- 
holders, an aperiodic roughness was obtained that was superimposed on the 
high-frequency screech. The larger V-type flameholder burned out after 
20 minutes of operation. The flameholders giving greater area blockage 
permitted a wider fuel-air-ratio range of operation. All the flameholders 
investigated screeched throughout most of their operating range of fuel- 
air ratio. 
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Fleimeholder position . - No appreciable change in screech frequency 
or amplitude was noticed when the flameholder position was altered to 
produce combustor lengths from 34 to 45 inches. At one position of the 
flameholder (giving a combustor length of 42 in.), however, a longitu- 
dinal mode of oscillation having a frequency of approximately 470 cycles 
per second appeared. This longitudinal oscillation did not replace the 
screech but merely accompanied and modulated it as shown in figure 17. 


Speculation Regarding Screech Mechanism 

The data thus far indicate that the gas in the burner does, in fact, 
execute a transverse oscillation during screech. Little has been said, 
however, of the means whereby the oscillation is maintained. Some source 
of energy is obviously reqiilred. Two sources of energy are available in 
the combustor: (l) the kinetic energy in the flowing gas stream, and 

(2) the chemical energy released during the combustion process. Of these 
two sources, the chemical energy is far greater. In order that this 
energy be made available to overcome damping and maintain the oscilla- 
tions at a high amplitude, it is necessary for it to be involved in some 
time- varying work cycle. The criterion proposed by Rayleigh (ref. ll) 
is that the heat input should occur at locations where the pressure varies 
and at a time when the pressiore is near its maximum. For a case where 
the oscillations are driven by the kinetic energy of the flowing stream, 
the heat release may play a part in intermittently storing and releasing 
this kinetic energy in such a manner that the oscillation is maintained 
(ref. 12). In either event, it is obvious that the oscillation may be 
driven only when the heat release undergoes a variation with time. This 
may be brought about in many ways. 

A coupling between the chemical reaction and the acoustic oscilla- 
tions could result from the effect of pressure and temperature on the 
chemical kinetics. The rate of energy release per unit volume in high- 
speed combustion equipment has been shown to increase with increase in 
pressure and temperature (ref. 13) . A relation of this sort can be used 
to predict the variation in heat-release rate resiilting from an adiabatic 
compression of the burning gases. A similar result obtains if the com- 
bustion zone is assumed to comprise many laminar flames again subjected 
to adiabatic compression. Both approaches Indicate an amount of energy 
available for driving the oscillation that is proportional to the mean 
local heat-release rate and to the square of the dimensionless perturba- 
tion pres Slice . 

Such an approximation was applied to the experimental burner re- 
ported herein. At typical conditions, perturbation pressure of 0.5, the 
available energy was found to be about 3 percent of the total energy in 
the fuel supplied. 
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The oscillation is damped, its energy removed, in at least four 
ways: (l) sound propagation upstream, (2) mass transport downstream, 

(3) absolution and scattering at the outer wall^ and (4) nonlinear damp- 
ing of the waves in the manner discussed in appendix B. 

Damping of types (l) to (3) depends on amplitude squared. Type (4) 
damping is a function of amplitude cubed. The rate of energy loss due 
to (l) and (2) for the experimental biarner was estimated to be from one- 
third to one-half the rate at which energy would be supplied. Damping 
of type (3) becomes of comparable magnitude only when the walls are made 
highly dissipative. The nonlinear attenuation (type (4)) is significant 
only when the amplitude is large. 

The foregoing crude picture enables the following predictions to be 
made as to the variation of screech tendency with operational parameters. 
The dimensionless amplitude p wovdd be an index of the difference be- 
tween the driving power and the "linear" damping terms h) to (3); hence, 
it would serve as an index of screech tendency. This tendency will in- 
crease with increases in mass flow, inlet temperature, flameholder block- 
age, and mean heat-release rate. The tendency will increase with pressure 
only if the local heat-release rate per potmd per second also rises. If 
the dimensionless amplitude p is constant, the amplitude P - will 

increase linearly with pressiire. These general predictions seem to be in 
agreement with the experimental findings . 

The previous concept of damping and driving assimies the heat release 
to be in phase with the pressure. If this is not true, the driving force 
would be reduced. Also neglected were the possible roles of flame wrin- 
kling and detonation. A flame-front disturbance, possibly cavised by the 
scattering of waves at the flameholder, certainly exists, as evidenced by 
the excellent schlieren photographs of reference 3 and by the data shown 
in figure 25. The res\iltant time-varying heat release may well be an 
Important constituent of the driving force. 

It is quite possible that in the environment of the screeching com- 
bustor, the heat-release rate may at times more nearly approach that of 
detonation rather t han ttiat associated with steady conflagration. This 
again coiold alter the energy driving the oscillation. Neither of the 
last-mentioned factors, however, would alter drastically the trends 
predicted by the use of the simplified picture. 

There are certain beneficial aspects of screech. One of these is 
the increased combustion efficiency accompanying screech. Another is the 
possible gain in the over-all heat cycle for a given thermal efficiency. 
This can occur because most of the heat is released at a higher level of 
pressure with screech. It is uncertain how much of this excess energy 
is available above that required to overcome damping of the oscillation. 

To retain the beneficial aspects of screech while minimizing the destruc- 
tive effects would demand a thorough and quantitative evaluation of con- 
trol techniques. 
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SUMMARY OF RESUTCTS 

The results of an investigation to develop instrumentation and study 
screech in a 6-inch-diameter simulated afterburner can be summarized as 
follows ; 


1. Probe microphones were developed that can be utilized to measxire 
the frequency^ relative amplitude^ and phase relations of pressure 
oscillations at various positions within a screeching combustor. 

2. Calibration data were obtained that make possible the determina- 
tion of the absolute value of pressure amplitude from probe-microphone 
readings. In obtaining these calibration data, a new theory was proposed 
to accovmt for the nonlinear atteniiation of high-ampllt\ide sound in 
tubes. . 

3. From measvirements of pressure-amplitude distribution, frequency, 
and relative phasing during screeching operation of a 6-inch-diameter 
simulated afterburner, the acoustic oscillations accompanying screech 
were found to consist of the first transverse (sloshing) mode in the hot 
gases downstream of the flameholder. 

4. The small-scale afterburner investigated screeched throughout 
most of its operating range. Screech was absent only at conditions 
approaching flame blow- out; screech disappeared just before blow-out as 
pressure was decreased or as air flow was increased. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 1, 1953 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report; 
a mean soimd velocity relative to the gas^ ft/sec 
c local sound velocity relative to the gas, ft/sec 

d combustor dieuneter, ft 
f frequency, sec"^ 

n Integer denoting particular- harmonic 

P local pressure, Ib/sq in. abs 

P^y mean pressure, Ib/sq in. abs 

p perturbation pressure, (p/Pg^y) - 1 

t time required for a wave to become saw-toothed, sec 

u local gas velocity relative to fixed coordinates, ft/sec 

V propagation velocity of saw-tooth shock front relative to fixed 

coordinates, ft/sec 

V wave -propagation velocity, relative to fixed coordinates, ft/sec 

X distance, ft 

X distance reqiiired for a wave to become saw-toothed, ft 

a dam p ing coefficient, ft~^ 

ttg effective over-aill damping coefficient, + Op, ft"^ 

cxjj. viscous damping coefficient, ft"^ 

Op saw-tooth damping coefficient, ft”^ 

a' d amp ing coefficient for wave of unit frequency, a' = o,/f^, ft~^ 

3 constant 

X wavelength, c/f, ft 
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Subscripts: 

e conditions obtained at probe entrance 
max conditions obtained at wave maximum pressure 

min conditions obtained at wave minimum pressure 

n integer denoting particular hannonic 

r conditions obtained at sound source 
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APPENDIX B 

DISTORTION AND DIMINUTION OF SOUND WAVES OF FINITE 
AMPLITUDE IN TUBES 

It has been established (ref. 9) that finite-amplitude sound waves 
tend toward a saw-tooth wave form as they propagate. A portion of a 
single saw-tooth wave has been examined in reference 14, which showed 
that, up to pressure ratios of 2.5, the sound speed computed by the 
isentropic relation 


1 

I = (1 + P)^ (Bl) 

differs by less than 1 percent from that determined from the Rankine- 
Hugoniot relation 


(B2) 


(B3) 

was found to vary less than 1 percent for shocks up to P/Pa,v “ 

(i.e. , p = 1.5) . 

Since Q is constant for rarefaction waves and approximately con- 
stant for shock waves, equations (Bl) and (B3) define the relation 
between u, c, and p. 

From the foregoing considerations, a number of useful relations 
can be derived for the propagation of sound in tubes. 

End effect . - At high sound levels, the variation in static pres- 
sure at the entrance of the tube will not be equal to the pressure var- 
iation in the medium from which the lng)ulse is taken. This results from 
the assumptions that (l) wave energy is conserved and (2) the incident 
sound has no particle velocity component in the direction of the tube. 
With these assumptions the energy equation reads 

2 ^2 ^2 
u + 5c = 5c 


1 
S' 

(when X = 1.4) . Also the Riemann invariant 

Q = 5c - u = constant 


c Rl + p) (7 + pH 

a L 1 + 6(1 +“p)J 


(B4) 
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Within the tuhe^ the speed of sound and the particle velocity are 
related hy the equation 


5Cg - Ug = 5a 

The speed of sound may he related to sound pressure hy 


(B3) 




+ p) = 1 




(Bl) 


Comhining equations (Bl) , (B3) , and (B4) yields the following expression 
relating the sound pressure in the prohe entrance with that in the 
source: 


Pr " Be ~ 14 ^e 


( 2 ) 


Transition to saw-tooth. - A harmonic distortion of finite sound 
waves results from the variation in sound speed between the compressed 
the rarefied portions of the wave, as indicated hy equation (Bl). 
Thus, any finite wave, regardless of its initial form, will approach more 
more closely the saw-toothed shape as time passes; that is, the 
leading portion of the wave will become a shock wave, while the rarefac- 
tion will become more and more gradual. The rate at which this transi- 
tion occurs will, of course, depend on the magnitude of the wave. 


If a train of sine waves he engendered in an initially stationary 
medium, the waves will become saw-tooth in form in a time 


t 


4(v - a) 

' max 


(B5) 


Since 


V = c + u 
max max max 


(B6) 


equations (Bl) , (B3) , (B5) , and (B6) may he combined to give 


* 

t = 


7A 


24ap 


(B7) 


max 


The distance the wave travels during transition is then 

7A 


X 


24p 


(B8) 


max 
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These relations should apply to waves propagating in ideal gases. 
If the waves are confined in a tube^ the factor of viscous damping must 
be considered. The ideal saw-tooth wave may be represented by the 
Fourier expansion 


2p 


P = 


max 


^ n 


sin 2nfnl 


(I - *) 


(B9) 


Since the viscous damping increases with increasing frequency, a physical 
wave may more truly be represented by terminating the series with the 
n^^ harmonic. The finite series will then exhibit a pressure maximum at 
a point 


- I iTTi 

if ft in equation (B9) is set equal to an integer. 

If the dan 5 >ing of the n^^ harmonic of a wave of frequency f fol- 
lows the relation 


-a 


'f“n^ 


P = P e 
•^n e 


the rate of decay of p^ with distance is 

'li’n , 

ar - - Pn"- 


m 


(Bll) 


(B12) 


Assuming a uniform rate of growth of this harmonic, equations (B8) and 
(BIO) yield 



P P 
n max 

7X n-1 

24 n+1 


(B13) 


The value of p below which p will no longer grow is found by 
max n 

equating the sum of equations (B12) and (B13) to zero, yielding 

~ 0.292an“(n-l)c ^ 3 ^ 

^max (n+l)f 

In free air m = 2 and a is small for the frequencies of interest. In 
fairly small tubes, m is l/2j when a s a' ,yF, equation (3) becomes 


Bmax 


0.292ca' A/^(n-l) 

\^(n+l) 


( 4 ) 
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If the viscous (Kirchoff) damping coefficient is applied to the 
wave during the transition from sine to saw-tooth hy the relation 


p = p e 
max e ^ max 


-dj^x 


(B14) 


the transition time and distance can he recalculated hy the use of 


. xn-1 P- 

+-^n+l " Jq 


dt 


(B15) 


This becomes, hy equations (Bl) , (B3) , (B6), (B14), and (B15), 

^ 24Pe^j]j3jj 3 , (n+l) 


t* = 


In 


n o^a -7Ao^(n-l) + 


(B16) 


or 


n 


= at 


n 


= ^ In 


24p. 


'e,max 


(n+l) 


24Pe,max(’^+l) ‘ 


or, for large n, 



e^max 


Be, max “ 0.292Xaj^ 


(B17) 


(5) 


Nonlinear attenuation. - If the original pressure wave is of suf- 
ficient amplitude to generate a large number of harmonics, a nonlinear 
attenuation results. This occurs as a result of the overriding of the 
shock front hy the following rarefaction impulse. In considering this 
attenuation it is assumed that the wave will remain approximately saw- 
tooth in form and that the mean pressure will he found halfway between 
the TngYimnTn an6 the minimum pressures. The Rankine-Hugoniot equation 


u - u . 

max min 


min 


rr/V-u . 

5 1 min 

6 \ c . 

\ min 


min 


V-u . 
min/ 


(B18) 


may he written in quadratic form and solved as follows: 

V = 0.6 u ^„y + 0.4 i Vi-i n + /\/o.36(u^g^ - ^^min^ *^in 


(B19) 


Equations (Bl) and (B3) may he combined to give 
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r 1 

1 7 

u = 5a [(l + p) 



(B20) 


which, with equation (B19) and the assumption of equal pressure and 
rarefaction magnitudes, gives 


V ~ ^max / , 15 \ /T ^ 31 2 

Thus, it is seen that little error would be introduced by assuming that 
the propagation velocity of the shock front is equal to the speed of 
sound in the undisturbed medium. 


The rate at which the saw-tooth wave diminishes because of the over- 
riding of the rarefaction upon the shock front can be expressed in terms 
of a nonlinear damping coefficient Op, defined as 


1 dp. 




max 

p to 
max 


dp. 


ap 


max 


max 

dt 


( 6 ) 


The equivalent expression, in terms of peak and average wave velocities, 
is 


a = — ^ (v - a) (B21) 

p aA max ' 

which may be combined with equations (Bl), (B3), and (B6) to give 




12f 

a 


(1 + 


p ) 

max 


1 


(7) 


This nonlinear damping is assumed to be accompanied by viscous damping; 
the effective damping coefficient for saw-tooth waves in a tube is taken 
as the sum of the two effects: 


1 

1 0f 7 “I Of 

ae - Cj. — (1 + - — (B22) 

The variation of peak pressure with distance can then be obtained by 

integration, again using the binomial expansion of (l + P mav ) 

discarding terms of higher order than p^ . This yields the expression 

Tno.x 
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-dp. 


cx^ = 


max ^ 12 f 


■msLX 


dx 


^ « P 


max 


3 

7 


p2 

■t'max 


) 


+ ai, 


(B23) 


from which 



(B25) 


where Cq is the constant of integration. For a given frequency and 
tube size, this equation may be plotted as in figiire 9. 


To summarize, there are three regions with distinct characteristics 
through which the sound wave travels in a probe-type microphone; the 
entrance, the transition region, and the saw-tooth region. Considering 
these regions separately, curves may be drawn for a given probe diameter 
that will represent the internal sound level as a f\mction of source 
level anri tube length. The sound level at the tube entrance is less 
than that of the source by the amount indicated from equation (2) . Dur- 
ing the transition (the length of this region is determined from eq. 
(5)), the linear (Kirchoff) damping coefficient applies (eq. (B14)). 

When the transition is complete, the damping coefficient of equation 
(B23) must be used. Such a set of curves is shown in figure 10. 
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TABLE I. - OBSERVED FREQUENCIES FOR SIMULTANEOUS VARIATIONS 
OF MASS FLOW, EQUIVALENCE RATIO, AND BURNER PRESSURE 
[^Inlet-gas temperature, 1660° rJ 


Frequency^ 

cps 

Mass flow^ 
lb/ sec 

Equivalence 

ratio 

Mean 

pressure, 

^av-’ 

Ib/sq in. 

Inlet 
velocity^ 
ft/ sec 

3350 

2.005 

0.9339 

19.78 

310.9 

3500 

2.080 

1.104 

20.78 

307.1 

3500 

2.080 

1.146 

20.78 

307.1 

3400 

1.990 

.9968 

20.28 

301.0 

3400 

2.452 

1.007 

22.28 

337.6 

3400 

2.787 

.9964 

24.28 

352.1 

3400 

3.166 

.9881 

26.28 

369.6 

3350 

3.474 

.9595 

28.28 

376.8 
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Figure 1. - Diagranmatic sketch of 6- inch afterburner installation. 
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Figure 2. - Diagranmiatic sketch of research combustor. 
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Figiire 3. - Diagrammatic sketch of air-atomizing spray bar. 
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Figure 5. 
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Sketch of ’’phase" probe (not scaled). 
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(a) Resonator for infinitesimals. 


Figvire 6. - Apparatus for measuring attenuation constants. 
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(b) Resonator for high amplitudes. 

Figure 6. - Concluded. Apparatus for measuring attenuation constants, 
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Figure 7 . - Attenuation constant a 
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Amplitude, p 

Figure 0. - Amplitude- dependent sav- tooth damping attenuation coefficient. 
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Mean-to-peak amplitude 
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Distance, x, ft 

Figure 9. - Integrated theoretical attenuation curve for nonlinear and Kirchoff damping 
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Figure 11. - Wave form of "infinite" probe microphone and reference microphone for 1080 cycles 
per seconds. Pressure, 0.40 atmosphere. 
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(a) Frequency, 1080 cycles per second. 
Figure 12. - Probe -microphone calibration. 
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Reference amplitude, P^ax 
(b) Frequency, 2400 cycles per second. 

Figure 12. - Continued. Probe-microphone calibration. 
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Figure 12. - Concluded. Probe-microphone calibration. 
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Figure 14. - Ionization gap. 
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Coolln@-vater inlet 


Figure 15. - Photocell prohe. 
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Figure 16. - Amplitude profiles across combustor for seven flameholder positions. 
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(a) Screech plus longitudinal resonance. 
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(h) Simultaneous trace of filtered major components in 
figure 17 (a) showing modulation of screech by low- 
freguency note. 


Figure 17. - Wave forms of screech and longitudinal resonance. 
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(a) Predominantly low-frequency spectrum; measured 4 inches downstream. 
Figure 18. - Examples of spectra at wall downstream of f lameholders . 
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(b) Predominantly high-frequency spectrum; measured 3/4 inch downstream. 
Figure 18. - Concluded. Examples of spectra at wall downstream of flameholders . 
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(b) First transverse mode. 



Figure 19. - Pressure contours and particle motions for fundamental 
modes of cylindrical duct. 
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Figure 21 . - Amplitude and phase distrihutions for two flameholder positions. 
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Figure 22. - Microphone voltages about circumference of screeching 
downstream of V-gutter. 
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Figure 23. - Pressure signal from "two phase probe microphones showing phase shift 
along a diameter. 
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Figiire 24. - Phase of pressure signal relative to signal at top 
of combustor 3/4 inch downstream of flameholder for three 
flameholder positions. 
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(a) Ionization, traces at top and bottom of combustor 3/4 inch from wall. 



(b) Ionization and pressure traces at top and 60° from top of combustor. 


Figure 25. - Ionization and pressure traces 3/4 inch downstream of flameholder. 
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Figure 26. - Reconstruction of pressure and flame-displacement 
sequence . 
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Figure 27 . - Peak-to-peak screech amplitudes measured 5^ inches upstream 
of flameholder for four different flow rates. 
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Flgxore 28. - Lower pressure limits for screech and burning. Equivalence 
ratio, 1. 
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